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Abstract 

We calculate the free energies F for ?7(1) gauge theories on the d dimensional sphere of radius R. 
For the theory with free Maxwell action we find the exact result as a function of d; it contains the 
term log R consistent with the lack of conformal invariance in dimensions other than 4. When 
the U{\) gauge theory is coupled to a sufficient number Aj of massless 4-component fermions, it 
acquires an interacting conformal phase, which in d < 4 describes the long distance behavior of 
the model. The conformal phase can be studied using large Nj methods. Generalizing the d = 3 
calculation in arXiv:1112.5342, we compute its sphere free energy as a function of d, ignoring the 
terms of order 1/Aj and higher. For finite Nf, following arXiv;1409.1937 and arXiv; 1507.01960, we 
develop the 4 —e expansion for the sphere free energy of conformal QED^. Its extrapolation to d = 3 
shows very good agreement with the large Nf approximation for Nf > 3. For Aj at or below some 
critical value Acnt, the SU{2Nf) symmetric conformal phase of QED 3 is expected to disappear or 
become unstable. By using the A-theorem and comparing the sphere free energies in the conformal 
and broken symmetry phases, we show that Acrit < 4. As another application of our results, we 
calculate the one loop beta function in conformal QEDg, where the gauge field has a 4-derivative 
kinetic term. We show that this theory coupled to Nf massless fermions is asymptotically free. 
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1 Introduction and Summary 

The four-dimensional Quantum Electrodynamics coupled to Nf Dirac fermions is an original model 
of Quantum Field Theory; its predictions have been verified experimentally with high accuracy. 
If the fermions are massless, then the theory is conformally invariant for zero charge e, but the 
interaction effects are well known to break the conformal invariance. They produce a positive /3 
function for e, which means that the theory becomes free at long distances. 

The physics of QED is different in d 7 ^ 4. Then the free Maxwell action is not 

conformally invariant [^, but the one loop fermion vacuum polarization diagram induces a scale 
invariant quadratic term proportional to 

, ( 1 . 1 ) 


which is in general non-local. For d < 4 this term dominates at long distances, and well-known 
examples of such “induced QED” are the Schwinger model in d = 2 and the conformal phase 
of QED 3 [^[^. In d = 4 — e the conformal QED^ theory may be studied using the e expansion, 
because the f3 function 
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3 (dvr)- 


+ 0{e^ 


has a weakly coupled IR fixed point at = deir'^/Nf 
operator dimensions in QED^^ was introduced in [m 


( 1 . 2 ) 


The e expansion of various 
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For d> A the induced term (1.1) becomes important at short distances (for d = 6 , 8 ,... it is a 
local higher derivative term). A perturbatively renormalizable theory in d = 6 based on the local 

^ In section 


4-derivative action was stndied 
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and, more recently, in 


12 


3.1 we will calculate the one loop beta function for this 6 -d theory coupled to Nf Weyl fermions; 


it has the negative sign which produces the asymptotic freedom. 

Among the important physical applications of the conformal QED is the theory in d = 3 
coupled to massless Dirac fermions and/or complex scalars. An early motivation to study QED 3 
came from work on the high temperature behavior of four-dimensional gauge theory [^. More 
recently, its various applications to condensed matter physics have been explored as well (see, for 



Specifically, we will work with the 17(1) gauge theory coupled to Nf massless 4-component 

Dirac fermions ipF The lagrangian of this theory has SU{2Nf) global symmetry, which is often 

2 

referred to as the “chiral symmetry.” In QED 3 the fine structure constant a = ^ has dimension of 
mass; this makes the theory super-renormalizable. At short distances we find a weakly interacting 
theory of massless fermions and photons, where the field strength Ff^i, has scaling dimension 3/2. 
The short distance limit of QED 3 is scale invariant, but not conformal. This is because the free 
Maxwell action \F^yF^^^ is not conformally invariant in three dimensions The lack of conformal 
invariance of the free Maxwell theory translates into the fact that its three-sphere free energy F 
depends logarithmically on the sphere radius R [28| . In section we will generalize this result to 
free Maxwell theory on and show that its free energy contains the term log R. We will refer 
to the short distance limit of QED 3 as the UV theory. The fact that its F value, F/jV) diverges is 
important for consistency of the RG flows with the T-theorem. 

As QED 3 flows to longer distances, the effective interaction strength grows and various inter¬ 
esting phenomena become possible. The one loop fermion vacuum polarization diagram induces a 
non-local quadratic term (1.1) for A^, which dominates in the IR over the Maxwell term [^. Due 
to this effect, the theory flows to an interacting conformal field theory in the large Nf limit where 
e^Nf is held fixed. In the GET the scaling dimension of F^i, is 2. The scaling dimensions of other 
operators can be calculated as series in t/Nf (see, for example, p9| - [^ ). 

A different possibility is the spontaneous breaking of the SU{2Nf) global symmetry dne to 
generation of vacnum expectation value of the operator (h is written using the 4-d 

or directly 


^Its Weyl anomaly coefficient a = —55/84 may be calculated either using a dual AdS approach 
in d = 6 


15 


13 14 
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notation for spinors "0* and gamma-matrices). This operator preserves the 3-d parity and time 
reversal symmetries, but it breaks the global symmetry to SU{Nf) x SU{Nf) x 17(1). This mecha¬ 
nism was proposed in |19j , where it was argued using Schwinger-Dyson equations to be possible for 
any Nj] however, for large Nf the scale of the VEV becomes exponentially small compared to a. 
Subsequently, modified treatments of the Schwinger-Dyson equations suggested that the chiral 
symmetry breaking is possible only for Nf < A^crit- The estimates of A^crit typically range between 
2 and 10 [7|[3§{^. 

It is widely believed that the QED 3 must be in the conformal phase for Nj > A^crit; but a nearly 
marginal operator may appear in the spectrum of the GET as Nf is reduced towards A^crit- This 
operator must respect the SU{2Nf) and parity symmetries of the theory, and natural candidates 
are the operators quartic in the fermion field^ (see also 


3^ ) 0 When the quartic operator is 
slightly irrelevant, it should give rise to a nearby UV fixed point; there is a standard argument for 
this using conformal perturbation theory, which we present in section We will call this additional 
fixed point QEDg. Eor Nf = A^crit it merges with QED 3 , and for Nf < Ncrit both fixed points may 


31 


become complex |32, 34,36, 371. In this “merger and annihilation of fixed points” scenario, for 


Nf < A^crit the UV theory flows directly to the broken symmetry phase. Alternatively, both fixed 
points may stay real and go through each other. Then the QED 3 fixed point continues to exist 
even after the appearance of a relevant operator; this relevant operator may create flow from QED 3 
to the broken symmetry phase. If so, the edge of the conformal window may be associated with 
the dimension of some operator in QED 3 becoming so small that it violates the unitarity bound. 
This would be analogous to what happens at the lower edge of the conformal window in the M = 1 


supersymmetric gauge theory 38 


We will attempt to shed new light on the transition from the conformal to the symmetry 
breaking behavior by using the E-theorem and performing more precise calculations of F. Here 


F = — logZ^s is the 3-sphere free energy 22,^ or, equivalently, the long-range Entanglement 


Entropy across a circle 24,25 . The theorem states that for Renormalization Group (RG) flow from 
fixed point 1 to fixed point 2, Fi > F 2 . A proof of this inequality has been found using properties 
of the Renormalized Entanglement Entropy in relativistic field theories [39] (see also [40] ). 

In order to apply the E-theorem to RG flows among different phases of QED 3 , it is important 
to know their E-values. This is especially challenging for the interacting GET phase of the theory, 
this calculation was performed using the 1/Nf expansion with the result 


In 


28 


Tconf — Nf 


/ log(2) 


3C(3)\ 1 /7rVf\ 1 , 


(1.3) 


^We are grateful to Z. Komargodski for informing us of this possibility. _ 

^ In the compact theory, monopole operators may also become relevant as one lowers Nf 35 ; however, these 


operators transform in non-trivial representations of the SU{2Nf) flavor symmetry, and so they are not expected to 
be generated along the RG flow if the UV theory has exact SU{2Nf) symmetry. Monopoles may still condense, i.e. 
they may acquire expectation values in the spontaneously broken phase. 
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The first term on the RHS is the T-value of Nj free Dirac fermions, A^j-Ffree-ferm- Even though 
Econf — AfjFfree-ferm grows without bound for large Nf, the E-theorem inequality Eyv > -Fconf is 
satisfied. This is because Euv is infinite due to the divergent contribution of the free Maxwell 
theory. In section we review the large N description of conformal QED and generalize the result 
(1.3) by computing Econf as a function of d. 

Since we will be quite interested in E^onf for small Nf, in this paper we will apply a dif¬ 
ferent approximation method 41,42 . This method consists of developing the e expansion of 


F = — sm(7rd/2)Fgd for d = A — e. It relies on the perturbative renormalization of the field 
theory on the sphere S‘^~^ and requires inclusion of counter terms that involve the curvature ten¬ 
sor |43| - |47| . Applications of this method to the Wilson-Fisher 0{N) symmetric CFTs have produced 
high-quality estimates of Eo(n) in d = 3; they are found to be only 2 — 3% below the F values for 
the corresponding free UV fixed points of these theories [4H|42j . 

In this paper, we will perform a similar e expansion for F of the conformal QED, building 
on earlier work which developed the perturbative renormalization of QED on S^~ 
calculation is presented in section]^ and our main result is 


47 


51 


This 


~ AT 1 . Nd 

Fconi — -^/-^free-ferm ^ 


+ 


3 l 7 r 

90 


- 1.2597e - 0.6493e^ 0.8429e'^ 


0.4418e2 0.6203e3 


N 


f 


N 


f 


f) (^■^) 


Extrapolating this to d = 3 using Fade approximants produces results very close to the large Nf 
formula (1.3) already for > 3, see figure 1 ^ 

Applying the E-theorem, we find that RG flow from the conformal to the broken symmetry 
phase is impossible when Econf < Fsb- This puts an upper bound on the value A^rit where the 
conformal phase can become unstable 26 . Using our resummed e expansion results for Econfj we 


find that the value of Nf where E^onf = Esb rather robustly lies between 4 and 5, and our best 
estimate is Nf ~ 4.4. If we restrict to integer values of Nf, this means that for Nf > 5 the QED 3 
theory must be in the SU{2Nf) symmetric conformal phase. Therefore, our results give the upper 
bound A"crit < 4. The same upper bound is obtained if we use the large Nf approximation (1.3) to 
Econf) which was derived in 28 . The results obtained very recently using the e expansion of quartic 
operator dimensions , as well as computations in lattice gauge theory 52,53 , are consistent with 
our upper bound. 


4 



















2 Free energy of Maxwell theory on 

The action for Maxwell theory on a curved manifold is 

S = I = ^J d^x^A’^ + Ri^ + , (2.1) 

where we have used F^j^y = Vfj^A^ — V^Af^ and [V^, = R^A^. On a round S'^ of radius R, 

we have Rp = and so the action is 

S = A, . (2.2) 

The partition function is given by 

where G is the volume of the group of gauge transformations. One way to proceed is to split the 
gauge field into transverse and pure gauge parl|^ 

A^ = B^ + d^cP, = (2.4) 


Following 28 , we have 


DA = DBD{d(j)) = T>.BT)>Y^det'(-V2) 

d+1 

vol(G) = 27ryvol(5^ I D'^ , vol(5'') = = ^dR'^ > 


(2.5) 


where prime means that the constant mode is not included. Then, the partition function can be 
written as 


^ ^ VdeT(-V2) r ^ 1 _ VdeF(-V2) 

27rv^vol(5'^) J 27ry^vol(5'^) ^’ 

( 2 . 6 ) 

where the subscript ‘T’ indicates that the determinant is taken on the space of transverse vector 
fields. 

The eigenvalues of the sphere Laplacian —acting on a transverse vector and corresponding 

"^Equivalently, one can use Feynman gauge by adding a gauge fixing term Lax = This gauge is more 

convenient for perturbative calculations when interactions with matter fields are included, and we will use it in Section 

El 
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degeneracies are known to be (see e.g. 


^ R‘^ 

For a scalar field, one has 

,(0) 


9l = 


r(£ + 2)r(d-1) 


>1. (2.7) 


X^'^) = (0) ^ (2£ + d _ l)r(^ + d _10 

^ + 9 i r(£ + i)r(d) ’ - 


( 2 . 8 ) 


In the case of the scalar field, £ = 0 corresponds to the constant mode which is to be excluded in 
our case. Using these results, the free energy of Maxwell theory on S”^, ^Maxwell = — log Z, can be 
written as 

Fm.„„ = ) Ei'f' + + log(2.vW)). 


e=i 


£=1 


(2.9) 

In dimensional regularization, the following results for the sum over vector and scalar degeneracies 
hold 


OO OO CO 

Esf’-i. E 4 ”’ = o, ^ E^r^-i- (2.10) 

e=i e=o £=i 

These can be obtained for example by evaluating the sums for sufficiently negative d where they 
converge, and analytically continuing to positive values of d. Using these regularized identities, one 


can readily extract the radius dependence of the Maxwell free energy (2.9) to be 

1 


^Maxwell = - ^ log(e^i2"-'^) + (d) , 


( 2 . 11 ) 


where (d) is a radius independent function of d (with simple poles at even d). In particular, 
we see that ^Maxwell —^ +00 in the short distance limit for d < 4. The function (d) can be 


evaluated in continuous d by computing the non-trivial sums in (2.9), as we explain below. 


We first find it convenient to rewrite the free energy in the following way 


F 


Maxwell — ^vector 2F)fiin—sc T Treasure ; 


( 2 . 12 ) 


where we have defined 




vector — 


_ 1 (1) i„^^(^ + l)(^ + d- 2), 




^ OO 

+ )E 4 "’m 


£=1 


Fmim —sc — „ 


Frr 


27re2i?2 

1 1 // i-^2\ ^ V'' (0) + d — 1) , 

- logdet (-V^) = 1^ — ^2 -’ 

= \og{2'K\J vol(5‘^)). 


£=l 


£=l 


l{l + d-l) 


), 


(2.13) 
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The grouping of terms in (2.12) is essentially equivalent to doing the calculation in Feynman gauge, 
where one has an unconstrained vector and a complex minimally coupled scalar ghost. To proceed, 
we use the identity 

(2.14) 


df 

log{y)= I 


'0 


Then, using the dimensionally regularized identities (2.10), one can rewrite the vector contribution 
as 




1 


dt 


vector 




li=i 


- ^ log(27re2). 


(2.15) 


Note that the radius dependence in Fvector, and the terms proportional to e *, have dropped out 


due to (2.10). The sum over I can now be evaluated analytically, leading to elementary functions 
of e“*. To perform the t-integral, it is convenient to use the identity 


1 


1 


t 1 — e 


— t 


due 


—ut 


(2.16) 


This allows for an analytical evaluation of the t integral, and after some algebra and using gamma 
function identities such as r(x)r(l — x) = 7rcsc('7rx), we arrive at the result 


Fvector — / du 

Jo 


IT 


(m + 1 — 3d{l + tt) + 2u{u + 2)) sin(—(d — 2n)) 


T (d — 2 — tt) T (1 + u) 


2sin(f)r(d) 


d- 2 


- ^ log(27re^). 


(2.17) 


(d-2)2-u2 

To evaluate the ghost contribution by similar methods, we can introduce a small regulator 


to deal with the zero mode, so that we can extend the sum over all modes and make use of (2.10) 


.^min—sc — 


dt 




£=0 


(2.18) 


Performing first the sum over i, using (2.16) and evaluating the t-integral, we obtain, after sending 
d —)• 0 at the en43 


F- --Cdu 

Jo 


,, . /TT/, r (d — u) r (u) 

(d — 2u) sm( —(d — 2n)) 


2sin(f )r(d + l) 2u 


( 2 . 18 ) 


We can now put everything together in (2.12) and obtain the radius independent part of the 


®We use log((5) = - + log(l + 5). 
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Maxwell free energy (2.11). We find 


F, 


( 0 ) 


Max, 


(d) = ^ log (27r(d - 


■ / dufdiu), 

sin(^) Jo 


( 2 . 20 ) 


where the form of fd{u) can be read off from the above results, and it is equal to 


TT 


r{d- 2 -u)r{l + u) sin(^)((i-2) 


fd{u) = —(u + 1 — 3d(l + u) + 2u{u + 2)) sin( —(d — 2u)) 


TT T(d — u)T(u) sin(^) 

- (d - 2u) ..n(-(d - 


2r{d) 


+ 


{d - 2)2 - u2 

( 2 . 21 ) 


Here the first line comes from the vector contribution (2.17), and the second line from the ghost 


contribution (2.19). Note that the the UV divergences of the free energy are fully accounted for by 


the overall sine factor in front of the integral in (2.20). 
Equivalently, in terms of F we have 


nl 


^Maxwell = bin(^)log(e^F^ "*) - bin(^)log(27r(d-l)2Hd) + J dufd{u), 


( 2 . 22 ) 


which is a finite smooth function of continuous d. 

As a test of this result, we can check that in d = 4 F reproduces the known value of the 


conformal anomaly a-coefficient for the Maxwell theory. From (2.22), we obtain 


F 


d=4 

Maxwell 


TT 7 , , , , q n , vr 31 

= — / du(l — u){vJ — vJ — llu + 12) = — • TV 


12 


2 45 


(2.23) 


corresponding to the correct a anomaly coefficient, a = fg (we use units where a = ^ for a 4d 
conformal scalar field). 

In other even values of d, the Maxwell theory is not conformal and F cannot be interpreted as 
an anomaly coefficient. Nevertheless, F still yields the coefficient of the 1/e pole in dimensional 
regularization, which fixes the coefficient of the curvature counterterm in the renormalized free 


energy. From (2.22), we find for instance 
vr 1271 


pd=6 


Maxwell 


2 1890 


pd=8 


TT 4021 


Maxwell 


2 6300 


j^d=10 
^ Maxwell 


TT 456569 


2 748440 


(2.24) 


The d = 6 result agrees with the value obtained in Appendix of 12 . For other even d values, we 
have checked that our results are in agreement with the coefficient of the logarithmic divergence 


for a massless spin 1 field obtained by zeta function methods on Euclidean AdS 2 n 55 
As a further check, in d = 3 we obtain the result 


-<d=3 


1 


Maxwell 2 




—+ - — —(2u^ + 3u^ — 23u + 12) cot(7rM) 
1 — u 12 






















(2.25) 


in agreement with psl. In d = 3, the Maxwell theory is Hodge dual to a compact minimally coupled 


scalar field. Note that from (2.19) we can read off the i^-value for a (non-compact) minimal scalar 
in d = 3, with zero mode removed, to be 


^min-sc = ^ log(7r) + ^ + log(i?) . 


(2.26) 


This result agrees with the one obtained in 128, 56 , and after carefully relating the radius of the 


compact scalar to the electric charge e, one can verify equality of the partition functions under 
Hodge duality. 

In d = 5, we find for F = —F: 


^MaLell 2 ^°^( 327 r 4 i?) 


+ [ du 
0 


= bog(^) + ^ + ^ 

2 16ir2 16 it ‘ 


- 7 ^ - — (6rt® — 35u^ -|- 275u^ — 486tt + 240) cot(7ra) 

w u2 - 9 240 ^ ^ > 

(2.27) 


It would be interesting to reproduce this result from a massless 2-form B 2 on S^, which is related 
by Hodge duality to the Maxwell theory. 


3 Conformal QED at large N 

The action for Maxwell theory coupled to Nf massless charged fermions in flat space is (in Euclidean 
signature) 

1 


Nf 


s = 






(3.1) 


2=1 


Here the fermions V’* are assumed to be four-component complex spinors. These correspond to 
Nf usual Dirac fermions in d = 4, while in d = 3 they can be viewed as 2Nf 3d Dirac fermions. 
In particular, in d = 3 the model has the enhanced flavor symmetry SU{2Nf). We define the 
dimensional continuation of the theory by keeping the number of fermion components fixed. In 
other words, we take 7 ^ to be 4 x 4 matrices satisfying { 7 ^, 7 ^} = 26^'^ 1, with trl = 4. All vector 
indices are formally continued to d dimensions, i.e. = d, 7 ^ 7 /g = d • 1 , etc. 

One may develop the 1/N expansion of the theory by integrating out the fermions. This 
produces an effective action for the gauge field of the form 


Seff = / d'^x^F^^'^F., - 


4e2 


I d'^xdS (^-A^{x)A'^{y){J,{x)Uy))o + 0{A^)^ , (3.2) 
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where 


(3.3) 


Jfj. = 

is the conserved U{1) current. Using the fermion propagator 


(V'*(x)V’j(0)) = -5) 


.r(i) 






d'^P 


27r5 (x 2)5 (27r)'^ p2 

the current two-point function in the free fermion theory is found to be 


^tpx 


(3.4) 


{Mx)M 0))0 = 


Cj = iV/trl 


(^ii) 


I d 

\ 2712 


(3.5) 


In momentum space, one find^ 

{Mp)M-p))o = [ d^xe-^P^J^{x)MO))o = -Cj- 


23-d^d/2Y (2 _ 0 

m 


9fj,u 


PfiPp 

p2 


(p2).-l 


(3.6) 


Thus, when d < 4, one sees that the non-local kinetic term in (3.2) is dominant in the low momentum 
(IR) limit compared to the two-derivative Maxwell term. Hence, the latter can be dropped at low 
energies, and one may develop the 1/N expansion of the critical theory by using the induced 
quadratic term 


>Jcrit — „ 


d‘^p 

(27r)“ 


A>^{p)A'^{-p){J^{p)M-p))o + 0{1/Nf) 


(3.7) 


Note that this effective action is gauge invariant as it should, due to conservation of the current. 

To compute the sphere free energy, we need to conformally map to S'^ and choose an appropriate 
gauge fixing. As in the previous section, we may gauge fix by splitting A^ = -|- 9^0, where 

= 0. Then, following the same steps as in (2.3), (2.6), the sphere free energy is given by 

F = A^/Ffree-ferm(d) “h ^ log detr ^ log det'(-V^) log |^27ryvol(S^^ 

where = —{J^Ju) is the non-local induced kinetic term, and Tfree-ferm is the contribution of a 
free four-component Dirac fermion (4l] 


Tfree—ferm(d) — 


sin(f )r(l + d) Jo 


du cos I — j r 


7r?x\ /l-|-d-|-rt 


1 -|- d — 


(3.9) 


The ghost contribution was already computed in the previous section, and is given in (2.19). To 


gp,iy 2 

^More generally, for a spin 1 primary operator of dimension A, one has {Jfj,{x)J,y{ 0 )) = Cj— —E'-nd 


\JP)} Fj r(A+i) 


2A-d PilP, 
A-1 p2 


-) ipr 
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evaluate the contribution of the transeverse vector, we hrst conformally map the current two-point 
function to the sphere of radius R, on which we choose the conformally flat metric 


ds^ = 


AR'^dx^dx^ 

(l + x2)2 


(3.10) 


Introducing the vielbein e^{x) = the two-point function for a spin 1 primary operator of 

dimension A can be written as 


{J^{x)Jy{y)) = Cje^{x)e]^{y) 


gmn _ 

\x-y\'^ 


s{x,y) = 


2R\x - y\ 


s{x,yY^ 

where in our case A = d—1, corresponding to a conserved current. The spin 1 determinant in 


(l + x2)l/2(l + y2)l/2 ’ 

(3.11) 


may be computed by expanding in a basis of vector spherical harmonics 13,28,57 . Splitting the 


vector Afj, in transverse and longitudinal parts, the spin 1 and spin 0 eigenvalues of = —{J^Ju) 
turn out to be, in the case of general conformal dimension A (see Appendix [A]) : 


(i)_ 2'^-2^7r''/2(A - l)r (f - A) r(£ + A) 

— —^J - 


r(A + i) 


r (^ + d - A) i? 2 A-d ’ 


^(0) _ d - 1 - A^(1) 


(3.12) 


Ar = 


A - 1 


1 ’ 


with degeneracies given in (2.7) and (2.8). For A = d — 1 the longitudinal eigenvalues vanish as 


expected, due to gauge invariance. The spin 1 contribution in (3.8) is then 

r{2-|)r(f)" 


llogdetA^ 


= - log ( Nf- 

2 2'^-27r2+^r(d)i?'^-2^ 


{i),„„^r(£ + d-i) 


+dEs;’i‘>8 


1=1 


r(£ + i) 


(3.13) 


where we have used the dimensionally regularized identity (2.10) to extract the constant prefactor 
in the eigenvalues, and used Cj = ANr . From this expression, we immediately see that 

V 2713 J 

the free energy contains a term ^ \og{Nf), independently of dimension. This can be traced back to 


the trivial constant gauge transformations on the sphere, or equivalently to ghost zero modes 13 


Note also that the radius dependence cancels out against the ghost and measure contributions 


in (3.8), as expected by conformal invariance. The remaining non-trivial sum may be evaluated 


directly for instance by using the integral representation 

logF(2) = 




-t 


(3.14) 


and following similar steps as described in the previous section. A compact form of the final answer 


for the sum is suggested by the results of 13 , where a formula for the change in F due to a 
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deformation by the square of a spin s operator of dimension A was computed using higher spin 
fields in AdSd+i with non-standard boundary conditions. For spin 1 , that result implies: 


- OO 


1=1 


T{e + A) \ 1 


T{i + d-A) 

r^-i 


, (0), (d-l-A r(£ + A) 

+ ‘°g( A-i r« + d-A ) 


1=1 

. . 2 ^ 


sin (^) r {d) Jo 


duu{d^ — 4u^) sin(7rn)r — l-t-n FI- — 1 — n .(3.15) 


d 


Taking carefully the limit of A = d — 1, and using 41,^ (note that the sum starts from i = 0 
here): 


1 


E4”’i°8 


£=0 


r(i + A) 


1 


r-A-; 


T{£ + d-A)J sin(^)r(d + l) Jo 

we finally obtain the result 

(i)i fr{i + d-i)\ 1^ ^r(d-i)' 


d 


du u sin(7ru)r I - + u F I - — u 


d 


(3.16) 


e=i 


— du 


{d — 2)‘^{d — l)u (4 + — (d — 2)^u^) 


sm 


7r{d—2)u'\ p / (d—2)(1—tt) ^ / (d—2)(l-|-it) 


16sin(^)r(d+l) 


+ 


1 


2(1 -u) 
(3.17) 


We have explicitly verified that this agrees with a direct evaluation of the sum using (3.14). 

Putting everything together, the final result for the sphere free energy F, or equivalently for 
F = — sin(^)F, takes the fornj^ 


1 7T(] ( N f \ 1 

F = NfFfree-ferm{d) - - sin( —) log ( “— 7 ^ + Ao{d) + 0{ — ) 


sm( ^) / ' Nf 


(3.18) 


where 


Ao{d) =-sin(:^) 


+ )log 


5E4‘>iog 

i=i 

25-2d^3(^_ 2 ) 


T{i + d-l) 

r{£ + i) 


^ OO 

log (1(1+d-i)) 


£=1 


(3.19) 


and the sums can be given the integral representations in (2.19) and (3.17). The resulting Ao{d) is 


^Note that, due to the factor log {—Nf/ sm(^)), the free energy is real for 2 < d < 4, it has an imaginary part 
for 4 < d < 6, then it is real again for 6 < d < 8, etc. This is essentially due to the fact that the Maxwell term yields 
a contribution — | log(e^i?'*“‘*) to F, and at the RG fixed point e* is positive for 2 < d < 4, negative for 4 < d < 6, 
etc. 
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a smooth, finite function of d which is independent of R and Nf. In d = 3 it evaluates to 


Aoid = 3) = 


(3.20) 


and so we find agreement with (1.3). For comparison to the perturbative calculation in the e 


expansion given in the next section, it is also useful to expand (3.19) in d = 4 — e. We find 


SItt 

Ao{d = 4 - e) =- 0.905e - 0.649316^ + 0.374025e3 + 0{e^) 

90 


(3.21) 


The leading term correctly reproduces the a-anomaly coefficient of the d = 4 Maxwell field, as 
expected. In the next section we will reproduce the terms to order from a perturbative calculation 
on 

Let us also note that in d = 2 we find 


Ao(d — 2) — —— 

D 


(3.22) 


59 


one 


. This result is exact 


corresponding to a shift of the central charge by —1. This is as expected, since in d = 2 we get 
the Schwinger model coupled lo Nj massless fermions; via the non-abelian bosonization 
finds that at low energies it is a CFT with central charge c = ‘iNf — 1 pO 
(all the 1/A^/ corrections to F should vanish as d —>• 2), and we will make use of it in Section]^ to 
impose a boundary condition on the Fade extrapolations of our e expansion results. A plot of the 
function Ao(d) is given in Fig. 


^o(cl) 



Figure 1: Plot of the smooth function Ao(d) from eq. (3.18). It has values F = 
—55/84) in d = 6, F = 3l7r/90 (a = 31/45) in d = 4, and F = —vr/G (c = —1) in d 


—557r/168 (a = 

= 2 . 
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3.1 Comments on d > 4 


In d > 4, one still formally finds a conformal electrodynamics in the large momentum (UV) limit, 


see eq. (3.6), but the corresponding CFT’s are non-unitary. For instance, in d = 6 the induced 


kinetic term (3.6) corresponds to the conformal spin 1 gauge field with Lagrangian L ~ 


13,15,62 . The a-anomaly coefficient for this conformal field can be extracted from our general 


result (3.18) setting d = 6 , which yields 


{F - A^/^free-ferm)Id=6 = ^ I (213u^ + 6 u® - 630 m^ + - ISdu^ + 314u - 120 ) 

240 do 

TT 55 , , 

=- 7 • (3-23 

2 84 ^ ^ 

corresponding to a = —1| (in units where a = ^ for a 6 d conformal scalar). This agrees with 
the result for the a-anomaly of the 6 d conformal spin 1 field, which can be obtained from one-loop 


determinants in AdS^ with non-standard boundary conditions 13,14 , or by a direct computation 


15 


. Note that this is not equal to the coefficient of the logarithmic divergence for a ordinary 


Maxwell field in d = 6 , eq. (2.24). As recently observed in 12 , this conformal spin 1 field is part 


of a non-unitary N = (1,0) conformal multiplet including a Weyl fermion with 3-derivative kinetic 
term and 3 conformal scalar fields, whose total a anomaly coefficient is o = — Hg, which turns out 
to be the value assigned by definition in to a AA = ( 1 , 0 ) vector multiplet in d = 6 . 

For finite Nf, one approach to the conformal QED in d > 4 is to use the d = 4 -|- e expansion. 


From the one-loop beta function (1.2), one sees that there are UV fixed points at imaginary values 


of the coupling. The large Nf limit considerations discussed above strongly suggest that these UV 
fixed points have a UV completion in d = 6 — e as the IR fixed points of the higher derivative 
renormalizable gauge theory 


1 


5 = I d'^x{ ^Ff,y{-V^)F>^'^ - + iA^W 


(3.24) 


where V’* are Vj 6 d Weyl fermions. To get an anomaly free theory, we may add Nf Weyl fermions 
of the opposite chirality, so that the model includes iVj 6 d Dirac fermions. The one-loop beta 
function for this theory can be computed by evaluating the correction to the gauge field propagator 


due to the fermion loop, which is given by (3.6) for general d. Expanding in d = 6 — e, one finds a 


pole that fixes the charge renormalization, and for the theory with Nf Dirac fermions, we obtain 
the beta function 

Nf 


2^ 1207r3 


+ 0{e^). 


(3.25) 


Unlike the case of QED 4 , this theory is asymptotically free in d = 6 . It would be interesting to 
compute the beta function for the non-abelian version of this model. By analogy with d = 4, we 
expect that in this case the pure glue should give a positive contribution to the beta function. 
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while matter gives negative contributions. The pure glue theory may then have IR fixed points for 
positive that could provide a UV completion of the Yang-Mills theory in d = 4 + e. Also, in the 
theory with gauge group SU(Nc) and Nf massless fermions, one may contemplate the existence of 
a conformal window directly in d = 6. 


4 Sphere free energy in the e expansion 


The action for massless QED in d = 4 — e on a general curved Euclidean manifold is 

Nf 

s = [ + aoW^ + boE + c^V}/{d - l)^) , 

(4.1) 


where are iVj four-component Dirac fermions, and we have added a Eeynman gauge fixing term, 
which we find most convenient for the perturbative calculation below. Here is the curved 
space covariant derivative (when it acts on fermions, it includes the spin connection term as usual). 
Finally, TZ denotes the Ricci scalar, is the square of the Weyl tensor and E is the Euler density: 

E = - 4K,,„K'“' + ■R?. (4.2) 


The action includes all terms that are marginal in d = 4, and cq, oq, 6 o, cq are the corresponding bare 
coupling parameters. Renormalizability of the theory on an arbitrary manifold implies that the 
divergencies of the free energy can be removed by a suitable renormalization of the bare parameters 
which is independent of the background metric. The renormalization of the electric charge is fixed 
by the flat space theory and reads, in minimal subtraction scheme [^: 


eo = /U 2 e 


4Nf ^ /8Nj ^ e® ^ f WONf ^ 88Nj 2Nf{22Nf + 9)' 


+ Hv^ + - 


3e (47r)2 v Se^ 


(47r)4'^l 


27e3 9e2 


27e 


where e is the renormalized coupling, and the corresponding beta function ii 



4:Nf 4Nfe^ 2Nf{22Nf + 9) 2A'/(4iV/(154iV/2808C(3) - 855)-h 5589)e^ 


^ 2 3 (47r)2~^ (47r)^ 


(4vr) 


243(47r) 


(4.4) 


®The terms of order e® that we have omitted from (4.31 can be reconstructed from (4.4 1 if desired. 
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Then, one finds an IR stable perturbative fixed point at 


e* = TT 



1 -- 


9 3(44iV/ + 207) 2 (2464iV2 U928NfC{3) - 45756iV/ - 62937) 


16iV. 


-ed 


512iV2 


-e"+ 


24576iV3 


e=^+ 0 (e^)) 


(4.5) 


The first few terms in the renormalization of the curvature couplings have been obtained in 
for Nf = 1, and in 47 for the general case. In our conventions, they reac0 


, Nf + 2 7Nf e2 

ttQ — fj. ( a + „„ ,,—^ + „„ 7 :;—^ + ... 


20e(47r)2 72e (47r)4 

, _ lliVf + 62 Nf /2Ar2 {imf + 9)Nf\ 

^ ~ 360e(47r)2 (A-k)^ ^ (iP R)8e / (47r)8 


360e(47 

_ f Nf ' 


(4.6) 


and the corresponding beta functions for the renormalized parameters a,b,c are 


^ Nf + 2 7Nf e2 

^““"“+20(4;^ + ^(4)^ + 
llA^f + 62 Nf 

= e 6 - .0 + ^ 


/3c = ec- 


360(47r)^ 
4.Nj e 6 
9 (dvr) 


2 (47r)® 

+ ... . 


(16iV/ + 9)A^/ e® 
27 X^) 


+ 


(4.7) 


We are interested in computing the free energy of the theory on a round sphere 5''^ of radius 
R, for which one has = 9,E = d{d — l)(d — 2){d — 3)/R‘^,TZ = d{d — l)/7?^. In particular, 
the renormalization of the Weyl square coupling qq will not play any role in this calculation. 
After renormalization, the sphere free energy F(e, b, c, fxR) is a finite function for any value of the 
renormalized couplings e, b, c. By standard arguments, it satisfies the Callan-Symanzik equation 


("li++'’4+-*4) 

As explained in |42| , it follows that to obtain the radius independent free energy at the IR fixed 

point we should set not only /3e = 0, but also the curvature beta functions (db = l^c = 0. The 

®The term of order e®/e in 60 is scheme dependent in the sense that it depends on the definition of the Euler 
density E in d = 4 — e. Our conventions for E in ( |4.2| differ from by an overall d-dependent factor. One can 
verify that the free energy at the fixed point is not affected by this convention dependence. 
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corresponding fixed point values in d = 4 — e are given by e = e* in eq. (4.5), and 
l/llfV/ + 62 Nf et (16Nf + 9)Nf e® 


6 * = - 


c* — 


e V 360(47r 
14N^ e 6 


2 (dvr)® 


+ 


27 


(dvr)^ 




(4.9) 


.. 9 (dvr)® 

and the free energy at the fixed point is the radius independent quantity 

d^conf (^) — . 


(4.10) 


Note that at the free field level, the effect of (4.9) is simply to remove the coupling independent 


part of the curvature terms. The fixed point free energy has then a pole due to the free field 
determinants, but F = — sin(7rd/2)T is finite in the d —)• 4 limit and reproduces the a anomaly of 
the 4d theory |411|42|. 


To perform the explicit calculations, we find it convenient to follow [^,64-66 and describe the 
sphere by flat embedding coordinates a = 1, 2,..., d + 1 satisfying = B?. In this 

approach, one also introduces Dirac matrices Oq of dimension 2^^/^ satisfying the Clifford algebra 
in d + 1 dimensions {oa, ct^} = 2Sab, a,b = 1,..., d + 1. In this embedding formalism, the vertex in 


(4.1) is given by 


Tairj) = ieoQab{v)oib, Qabiv) = bab - 


riallb 

’ 


(4.11) 


where we have rescaled the gauge held to bring the coupling constant in the vertex. One advantage 
of using embedding coordinates is that the propagators take a relatively simple form |48| . The 
photon propagator in the Feynman gauge is 


Dabim-im) = babDim^m) = ^ab 
and the fermion propagator is 


r(d-2) 


(47r)2ii'^-2r(f) 




(4.12) 


= -S' 


r(f)a- (r?i -m) 


^ 27 r 5 \vi -m\‘^ 

Introducing the integrated n-point functions 

/ n 

k=i 


(4.13) 


(4.14) 
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the free energy is given by 


^QED^ = A^/-Ffree-ferm(d) + “ 2 "*) + 7r,^oG2 “ 7760^4 + • • • 


2 ! 


4! 


+ ndR^-^{d{d -l){d- 2){d - 3)bo + d^co) , 


(4.15) 


where = 27r t /r(^^) is the volume of the unit sphere, the free fermion free energy is given in 


(3.9), and we have separated out the coupling dependent part — | log(eQi?‘^ '^) of the free Maxwell 


free energy, see eq. (2.11), (2.20). This term plays an important role in the renormalization proce¬ 


dure upon using (4.3), and its presence is necessary for cancellation of poles and to obtain a radius 


independent free energy at the fixed point. 

The technical details of the calculation of G 2 and G 4 are given in Appendix To the order 
needed here, we find that their e expansion is given by 


G2 = 

G4 = 


„i ! ' 4(5 + 3(log(4ifl2)+7)) ISir^ + 124+ 4{5 + 3(log(47rR'^)+7))“ , .n, 2.4 

^^4^ +- liv -+- iw -' + ^h 


- 72A^j(5 + 3 (log( 47 ri? 2 ) +7)) + 4(77 + Dvr^) + 9iV/(72C(3) - 47)) + 0(e). (4.16) 


Plugging these results into (4.15), as well as the coupling renormalization (4.3) and (4.6), we 


hnd that all poles indeed cancel for arbitrary e, b, c. In particular, our calculation provides an 


independent check on the curvature counterterms (4.6) to order 


We can now compute the free energy at the IR fixed point by plugging in the critical couplings 

-P(;onf (e) = -^QED^ («*, b*, C*, ^R) 


(4.5), (4.9). Defining 


(4.17) 


we hnd 


^conf = A^/Tfree-ferm(d) + F^Lid) + ^ log 

+ (40 + 24(7 + log(4.)) + ^) ,(; + (>^^ + f - + 4 7 0(.»). (4,18) 

Note that the result is indeed independent of the radius R, consistently with conformal invariance 


and the Callan-Symanzik equation (4.8). 


( 2 . 20 ) 


Using the explicit e expansion of the free Maxwell contribution, which can be obtained from 

(4.19) 




we then obtain our hnal result for Fconf = — sin(^)Fconf given in eq. (1.4). As a test of this result. 
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one 


can verify that to order Nj it precisely agrees with the large Nj prediction, eq. (3.18) and 


(3.21). 


5 Fade approximation and the F-theorem 


A novel feature of the result (4.18) compared to the sphere free energy for the 0{N) Wilson-Fisher 
fixed points [^, is the appearance of the log(e) behavior in d = 4 — e. This makes it difficult 
to apply standard resummation techniques like Fade approximants. To circumvent this problem, 
we isolate the logarithmic term ^log{Nf/e) in F, which we treat exactly, and perform a Fade 
extrapolation on the function 


~ ~ 1 Trd Nf 

SFd{Nf) = Fconf - A^/Ffree-ferm + ^ sin( —) log(-^) 


(5.1) 


In d = 2, the IR fixed point of QED^ is the Schwinger model with 2Nf massless two-component 


Dirac fermions. In the infrared it is described, via the non-abelian bosonization 59 , by the level 1 
SU{2Nf) WZW model 60,61 . This is a CFT with c = 2Nf — 1 corresponding to F = f (2A^/ — 1). 


Therefore, it is natural to use a two-sided Fade approximant subject to the constraints: 


-f, d = 2, 

^ - 1.2597e - 0.6493e2 -h 0.84296^ -h ^ = 4 _ ^ _ 

(5.2) 

This allows us to use Fade approximants Fade[m,n] of total order 4. The results for d = 3 using 
these two-sided approximants Fade[ 2 , 2 ] Fade[i 3 ] given in table For comparison, we also 
present the results using one-sided approximant Fadep 2 ] obtained without assuming the boundary 
condition at d = 2 (we see, however, that its agreement with the large Nf expansion is not as good 
as that of both two-sided Fade approximants). We also plot the Fadeji 3 ] approximant for different 


Nf 

1 

2 

3 

4 

5 

6 

10 

50 

100 

Fade[ 2 , 2 i 

- 

-0.1512 

-0.1284 

-0.1237 

-0.1223 

-0.1218 

-0.1217 

-0.1230 

-0.1233 

Fade[i,3i 

-0.2743 

-0.1462 

-0.1284 

-0.1228 

-0.1204 

-0.1192 

-0.1176 

-0.1172 

-0.1172 

P Q^d-Gaverage 

- 

-0.1487 

-0.1284 

-0.1232 

-0.1213 

-0.1205 

-0.1196 

- 0.1201 

-0.1203 

Fade[i 2 i 

- 

-0.1856 

-0.1259 

-0.1072 

-0.0986 

-0.0937 

-0.0861 

-0.0799 

-0.0793 

e-expansion 

-0.7148 

-0.2113 

-0.1049 

-0.0632 

-0.0418 

-0.0291 

-0.0074 

0.0122 

0.0141 


Table 1: Various Fade approximations and the unresummed e-expansion of 6Fd{Nf) at d = 3. The 
two-sided approximants Fade[ 2 ^ 2 ] Fadep 3 ] are obtained assuming the value — | at d = 2, while 
Fadeji 2 ] does not use this assumption. Row 3 is the average of the two-sided approximants, i.e. of 
the first two rows. At large Nf we expect to find dFi =3 = 5 log(f) ~ —0.1208. 

Nf in figure]^ as a function of 2 < d < 4. In figure]^ we plot Tconf-QEDg — AijFfree-ferm comparing 


dFd{Nf) = < 


19 




























6~F,{Nf) 




Figure 3: Comparison of the Fade resummation of the e expansion, and the large N result (1.3) 
for the sphere free energy of conformal QED 3 . 


It is interesting to compare this result for the conformal phase of the theory with the F-value in 
the broken symmetry phase. The latter contains 2N'j Goldstone bosons and a free Maxwell field, 

At long distances, each of these fields is described by a conformally 


which is dual to a scalar. 


10 


coupled scalar field. Therefore, after the chiral symmetry breaking 

Fsb(JV/) = (2JV| + 1)(‘^-^), 


(5.3) 


To study if the F-theorem allows flow from the conformal phase to the phase with global 


In compact QED this scalar is compact 67 , and it develops a vacuum expectation value. If there are no monopole 


operators in the action, then the topological U{1)t symmetry is spontaneously broken and there is a massless scalar 
degree of freedom in the IR, just as in the non-compact case. We thank Z. Komargodski and S. Pufu for discussions 
on this. 
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symmetry breaking, we define the function 


mf) = FconiiNf) - ( 2 iV| + 1 )(^ - . (5.4) 


Its plot obtained using the Padeji 3 ] with d = 2 boundary condition is shown in Fig. One can 
also consider the corresponding function using the large Nf expression (1.3) for Pconfi this gives 
results that are close to those shown in Fig. 



Pade[i,3i 

Pade[2,2i 

Pade[i,2i 

Nf,c 

4.4204 

4.4180 

4.4530 


Table 2: Estimates of 
approximants. Padeji 3 ] and Pade 
only uses data from the d = 4 — e expansion. 


which is the solution of F^oni = Tsb, obtained from various Pade 
5 ( 2 , 2 ] use the d = 2 boundary condition in (5.2), while Pade[i 2 ] 


MNf) 



Figure 4: Plot of A(1V/) = Fconi{Nf) - FsB{Nf), using Pade[i_ 3 ]. 

The plot in Fig. [^implies that the RG flow from conformal to symmetry broken phase becomes 
impossible for Nf between 4 and 5. This value of Nf if found by solving the equation Fconf = Psb; 
it provides an upper bound on the integer value of Nf where the conformal window may become 
unstable: Acrit < 4. 

Now let us treat Nf as a continuous parameter and discuss the implications of conformal 
perturbation theory for the phase structure of the theory. As explained in the introduction, we 
expect that for Nf near Acritj a FU{2Nf) invariant quartic operator is nearly marginal, and we 
can work perturbatively in the small parameter 5 = A — 3. The beta function for A, the coefficient 
of the quartic operator, has the structure /3 a = dA + AA^ + O(A^). Thus, in addition to the QED 3 
fixed point at A = 0, we find a nearby fixed point at A* = —6/A. For Nf > Ncrit, this is a UV 
fixed point. It is another SU{2Nf) invariant CFT which we could call QED 3 . Its existence for Nf 
slightly above Ncnt is guaranteed by the conformal perturbation theory. It also exists for large Nf, 
where it is a double-trace deformation of QED 3 . Therefore, QED 3 may exist for all iVj > Ncrit- 

A commonly discussed scenario is that the QED 3 and QED 3 fixed points merge at Acrit and 
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(a). Nf > (b). Nf < 

Figure 5: Schematic picture of RG flows for Nf > A^crit (a) and Nf < A^crit (b). The QED 3 
and QEDg fixed points merge at A^crit and acquire small imaginary parts for Nf < A^crit- la the 
latter case, the interacting conformal behavior is no longer possible, but the RG flow from the 
UV can “hover” near the complex fixed points before running away to large quartic coupling and 
presumably leading to the broken symmetry phase. 


acquire small imaginary parts for Nf < A^crit [^[^[^[^. This means that the interacting 
conformal behavior is impossible for Nf < A^critj but the RG flow from the UV can “hover” near 
the complex fixed points before running away to large quartic coupling and presumably leading 
to the broken symmetry phase (see figure . During the hovering F can be made parametrically 
close to Tconf (Vcrit)- This is why the F-theorem requires Fconf(A^crit) > Tsb (a similar argument 
involving the continuity of F was given in [^). As we have seen, this gives a rather stringent 
bound Vcrit < 4.4. 

An alternate possibility is that both fixed points stay real and go through each other. Then the 
QED 3 fixed point continues to exist even after the appearance of a relevant operator; this relevant 
operator may create flow from QED 3 to the broken symmetry phase. Therefore, the F-theorem 
bound on Vcrit is the same as with the “merger and annihilation” scenario. 
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A Eigenvalues of the kernel 

On the S'^ in stereographical coordinates the kernel K^i, = —has the form 


K^u{x,y) = -Cj- 




(<5/.. - 2 


{x-y)^{x-y)t, 

\x-y\^ 


) 


(l+x2)(l+y2) 

We need to decompose the kernel on a sum of vector spherical harmonics: 


(A^l) 




ds) 


(A.2) 


^,ra s 


where s denotes different types of vector spherical harmonics, i is the principal angular quantum 
number and the range of m for given £ and s is (see (2.7)-(2.8)). The eigenvalue has 


no m 


dependence because of rotational invariance. Because of vector spherical harmonics are orthonormal 
to each other we havtF^ 




yd — id 


(«) 




(A.3) 


We first consider longitudinal vector harmonics 

■w(O) (^\ — 


MAmV ; + d-l) 


(A.4) 


where Ygmix) are usual scalar spherical harmonicsIntegrating by parts in (A.3) we get 




A Q/] m 


£{£ + d-x)j ^ 

where we sum over m and divide by the degeneracy We use that 


Dm ^ ' 


1 Cf-^^^^{1-S^{x,y)/2R^) 


9i m 


^(d-l)/ 2 (i) 


(A.6) 


^We assume that / d‘^x^/g^Y^y^{x)Ypl^i\x) = Su'Smm'Sss'- 


assume that V'^Yim = —i{i + d — 1)R "^Yem and f d‘^Xy/lRYemix)Y£j^,{x) = Su'Sr, 
^^We are allowed to do this sum because on the l.h.s doesn’t depend on m. 


^^Notice that the Y{m(x) harmonics have a factor of R which is consistent with the formula (A .61 because 
volfS^*) = 27r^i?‘^/r(^) and have a factor of R-'^^^+^. 
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where is the Gegenbauer polynomial. Now due to rotational invariance we may choose 

y = 0 and fincf^ 


V„V„A'»“'(x, 0) = - d + 1){2 - d + 2A + 

Therefore we obtain 

Cj{A-d + l) vol(5'^-i) 


(A.7) 


- - 


2^R^^-de{d + i-l) 


j dz{l 


+ {1 - z) {2 +2A-d-{1 - z){l +A-d))Cl‘" "^'\z), (A.8) 


(rf-l)/2. 


where we have changed the variable = (1 — 2 ;)/(l + z'). Calculating the integral we gel^® 

2'^-2^7r‘^/2(A-l)r(f-A)(i-A-l r(; + A) 


aS“> = -C,: 


1 


r(A + i) 


A-l r((i + Z - A) i?2A-d 


(A.9) 


The eigenvalues for transverse spherical harmonics can be easily found for A = d — 1, using 


the embedding formalism. In this case the kernel has the form 48 

R-^ 


Kabidl^m) = -Cj 


-Pr 


^cb 


(2d-4)(d-l) |??i - ’ 


(A.IO) 


where the operator Pac = \L'^5ac + LadLdc - {d- 1)Lac acts on r]i, and Lab = Va^j^ - Vb^ and 
L^ = LabLab ■ Now using the decomposition 48 


^cb 


\rii-y2\‘^d-i Z 


EE( 2 fi)' 


m s 


-d 2) (g) / NyW / X 

r(d - 2 )r(.^ + 2 ) b,lm\RV 


(A.ll) 


and the property of the operator Pab'- PabYb/m = 0 for s / 1 and PabYb^ln = 
we find 

Kab{vi,m) = '^Cj 

t,m 


( 2 fl)--.n<- + !)«■ + - 2)r(2 - j)r« + d - 2) ,1) ,1, 

■. R\2d - 4)(d - i)r(£; - 2)'r(< + 2)''' U.(mWiai,,,™W2), ia.i2) 


^®Notice that we can not just take j/ = 0 in (T.5I without having the sum over m. 

I6u u x-u • + 1 fi J n \ ( iNr2'’+"+5r(/J+i)r(i2+i)r(2i2+or{/9-i2+§) 

Here we used the integral J_^ dz(l + z) 2 (i _ z^Ci (z) = (-1) - ar( 2 i.)r(B-^-i+i)r(B+^+e+i) 

be found with the help of 
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r!r(2i2)r(/3-,2-r+|)r(/3+,2+r+|) 

7.311.3 and the relation C( (z) = (—1)^C/ {~z)- Also we used that ^ (1) = 
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therefore for A 


( 1 ) 


we get 


(1) ^ 7rir(2-|)r(£ + d-i) 1 

^ ^ 2'i-^r{d)T{£ + 1) ’ 


which coincides with (3.12) for A 


d-1. 


(A.13) 


B Calculation of G2 and G4 


Using the propagators and vertex in (4.13),(4.12), (4.11) we find for the two-point function 


G2 = -Nf I d'^riid‘^r]2tT{Qab{rji)<^bSir]i,'n2)Qcd{'n2)adS{r]2,'ni))Dacim,V2) ■ (B.l) 


After calculation we obtain 

(A^tri)r(|)r(d-2) f ^d-2 + ^im-v 2 ?) 


G2 = - 


2d+2^ 2 


d rjid ri2- 


\m - 


17 n .7 0^1 im - V 2 ?. 

2Fi{l,d-2,-,l - ^ -, 


(B.2) 


Due to rotational invariance we can put r/2 to the north pole of the sphere and get in stereographic 
coordinates 


(A/trl)ii^-‘^r(d - 2) 


G2 = - 


f 


dx 


d- 


)^U(M-2,^. ' ). (B.3) 


(167 r)Vr('^) Jo (l + x2)x‘^-U l + x2 

Now introducing the variable z = 1/(1 + we find 
{Nftrl)R^-<^T{d-2) ^ 


(1 -h x^)' 


d 


<^2 = ' / dzz^ C^-z) Gd-2z)2Fi{l,d-2,-,z). (B.4) 

2(167r) —r(^^) Jo 2 


The integral can be calculated exactly and we obtairp] 

{Nfill)R^-'^Tid - 1) 


G2 = - 


(B.5) 


(d-3)4^(47r)^sin(f )r(^^) ‘ 

After doing combinatorics for G4 we find that it consists of the sum of 34 integrals of the form 


14(01,..., oe) = 


r(|)N4 /i 22-dp('^_ 2)s 2 


2712/ V (47r)2r(f) 


n<i‘ 

2=1 




s(m> %)^“ 2 s ( 7 ^ 3 , 7/4)2035(7^^^ r74)2“4s(7^^^ 7y3)2a55^^2, r?4)2«6 ’ 

(B.6) 


^G2 was also computed in ISOlpl by different methods. 
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where we used the exact form of the photon propagator (4.12). The integral I 4 , is represented 
diagrammatically in figure 





2 

2 a 2 

3 


Figure 6 : The diagrammatic representation for the integral 14 ( 01 ,..., og). 

The next step is to write the integral in the Mellin-Barnes (MB) representation and then use 
the Mathematica program [^[70] to find I 4 as a series in e . For the photon propagator we use 
the Mellin-Barnes representatiorj^ 


= - 


sm 


(f) 


1 


r*+ 2 oo 


dzVi-zni + z)T{d - 2 + z)ril - ^ 


2'^7r 2 2 27ri 


(B.7) 


Then using methods similar as discussed in 142] one can write the general MB form for the integral 


14(01,..., Og) — 


sin2(^f)r(f)4 


-{2R) 


2d-H4-2ELl“i 


(2^ / fldziT{-zi)T{-Z2)T{zi + l)r{z2 + 1) 


2^+7^^ I (5d-|-3) Y ( 

T{d -2 + 2 ;i)F (1 - ^ - zi)T{d -2 + Z2)T{1 - ^ - Z 2 )T 2 {d - 0145 zi, oi, 05 - zi\z 3 , Z 4 , 2 : 5 ) 


Fo(d - 0235 + zi- Zi,d- oi3g -h 2:2 - Z3, 03 - 2:5), 


(B. 8 ) 


where amnk... ^ o-m + cin + o,k + ■■■ and F-blocks are 

^ r(| - 6 )F(oi + b - l)T{a2 + b - l)T{ai + 02 + 6 - d) 

Fo(ai,a 2 , 6 ) = ^^^- 


r(|) F(oi)r(a 2 )F(oi-h 02 + 26 - d) 

^d /2 r(-2i)F(o + bi+b2 - Ya=i Zi) 


F 2 (a, 61 , 62121 , 22 , 23 ) = 


r(o)F(6i)F(62)F(d-o-6i -62) 

X r(d — a — 26i — 262 — 2i — 22 — 223) . 


-r(6i + Z1+ 23)r(62 + 22 + 23) 

(B.9) 


For some values of ai,...,ag the MB representation (B. 8 ) is divergent or exactly zero for any d 


(due to the term r(0)). To handle this problem we used an additional regulator 6, so we were 
^®We use the formula 2^1(0, b, c; x) = r(a)r(6)r?c-a)r(c-i,) du I-i^ dzT{-z)T{a + z)r{b + z)T{c-a - b - z){l - x)^. 
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consistently calculating the integrals 


..., 06 ) = h{ai + 6, ...,06 + < 5 ). 


(B.IO) 


Each integral out of 34 in G 4 depends on 5 and e. But the G 4 itself is free of 6 and depends only 
on e as it should. 

In general MB approach gives a result in terms of a series in e and 5 and each coefficient of 
this series is a sum of convergent Melling-Barnes integrals, which are independent of e and S and 
can be calculated numerically. The numerical result is usually equal to some exact combinations 
of constants like 7r^,('(3), e.t.c. To get an exact answer we used integer relation search algorithm 
PSLQ . The final result reads 


G4 = 


Nf 


(8A^f(5 + 3(log(47ri?^) + 7)) — 18) 1 / 

+ ^ ^- - - - + (l6iV| (5 + 3(log(47ri?2) + j)y 


1227r'^e 


- 721V /(5 + 3 (log( 47 ri? 2 ) + ^)) + 4(77 + 97 r^}Nf + 9A^/(72C(3) - 47)) 


(B.ll) 


Note that in this calculation of G 4 we have used d = 4 — e and N = iV/trl = 4iV/. 
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